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ABSTRACT 

Using a sample of optically-selected quasars from the Sloan Digital Sky Sur- 
vey, we have determined the radio-loud fraction (RLF) of quasars as a func- 
tion of redshift and optical luminosity. The sample contains more than 30,000 
objects and spans a redshift range of < z ^ 5 and a luminosity range of 
—30 ^ Mj < —22. We use both the radio-to-optical flux ratio {R parameter) 
and radio luminosity to define radio-loud quasars. After breaking the correlation 
between redshift and luminosity due to the fiux-limited nature of the sample, 
we find that the RLF of quasars decreases with increasing redshift and decreas- 
ing luminosity. The relation can be described in the form of log(RLF/(l-RLF)) 
= bQ+bz log(l+2)+6M(^25oo+26), where M2500 is the absolute magnitude at rest- 
frame 2500 A, and b^, bM < 0. When using i? > 10 to define radio-loud quasars, 
we find that bo = -0.132±0.116, b, = -2.052 ±0.261, and 6m = -0.183±0.025. 
The RLF at z = 0.5 declines from 24.3% to 5.6% as luminosity decreases from 
M2500 = -26 to M2500 = -22, and the RLF at M2500 = -26 declines from 24.3% 
to 4.1% as redshift increases from 0.5 to 3, suggesting that the RLF is a strong 
function of both redshift and luminosity. We also examine the impact of fiux- 
related selection effects on the RLF determination using a series of tests, and 
find that the dependence of the RLF on redshift and luminosity is highly likely 
to be physical, and the selection effects we considered are not responsible for the 
dependence. 
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1. Introduction 



Although quasa rs were first discovered by their radio emission (e.g. lMatthews fc Sandage 



1963 



Schmidtlll963l ). it was soon found that the majority of quasars were radio-quiet (e.g. 



Sandage llQGsl ). Quasars are often classified into two broad categories, radio-loud and radio 
quiet, based on their radio properties. There is mounting evidence that the distribution of 
radio - to-optical flux ratio for optical l y-selected quasars is b ir nodal (e.g. 



Kellermann et al. 



19891 : Miller. Peacock, fc MeadI Il990l : IVisnovskv et al.l Il992: llvezic et al 



20021). although 



the existence of the bim odality has been questioned (e.g. iCirasuolo et al. boosl . but see 
also llvezic et al.l (l2004al ) for a response). Radi o-loud and radio-quiet quasars are pro bably 
powered by similar physical mechanisms (e.g. Barthel 1989 : Urry fc Padovanil Il995l ). and 



their radio properties are also correlated with host galaxy properties, central black hole 



masses, black ho le spins, and accretion rates (e.g. iBaum et al. 



Baum et al. 


1995; 


Urrv & Padovani 


1995 



Best et al. 



20051). Radio-loud quasars are likely to reside in more massive galaxies (e.g. 



Peacock. Miller, fc Longaii 



holes (e.g. 



Laoi 



200C 



1986 



Lacy et al. 



Bes 



2001 



et al. 



20051), and harbor more massive central black 



McLure fc Jarvisll2004l ). than do radio-quiet quasars. 



Roughly 10%— 20% of all quasars are radio-loud (e.g. iKellermann et al.lll989l : lUrry &: Padovani 



19951 : llvezic et al.ll2002l ). However, the radio-loud fraction (RLF) of quasars may depend on 
redshift and optical lumi nosity. Some studies have found that the RLF tends to drop with 



increasing redshift (e . g. iPeacock. Mi 



Visnovskv et al 



1992 



Goldschmidt et al. 



Schneider et al 



1999 



and luminosity (e.g. iHooper et al 



er. fc Longairl Il986l : iMiller. Peacock, fc MeadI ll99ol; 



Cirasuolo et a. 



1995 



19921) an d decreasing luminosity (e.g. |Padovani.l993, : 



2003r). or evolves non-m onotonically with redshift 



the RLF does not differ sig nificantly with reds hift (e.g. 



2000l : ICirasuolo et al.ll2003l ) or luminosity (e.g. iBischof fc Becker 



Bischof fc Becker 1997). while others showed that 



Goldschmidt et al. 



1997 



199 



i 



Stern et al. 



Stern et al 



2000! 



Fr om a sampl e of 44 72 quasars from the Sloan Digital Sky Survey (SDSS; lYork et al. 



2000|), llvezic et al.l (|2002l ) found that the RLF is independent of both redshift and optical 
luminosity when using marginal distributions of the whole sample; however, they noted 
that the approximate degeneracy between redshift and luminosity in the SDSS fiux-limited 
sample may cause individual trends in redshift and luminosity to appear to cancel. By 
stacking the images of the Faint Images of the Radio Sky at Twenty-cm survey (FIRST; 
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Becker. White. &: Helfandl Il995l ) , IWhite et al.l (120061 ) were able to probe the radio sky into 
nanoJansky regime. They found that the median radio loudness of SDSS-selected quasars 
is a declining function with optical luminosity. After correcting for this effect, they claimed 
that the median radio loudness is independent of redshift. In this paper we use a sample of 
more than 30,000 optically-selected quasars from the SDSS, and break the redshift-luminosity 
dependence to study the evolution of the RLF. We will find that there are indeed strong 
trends in redshift and luminosity, and that they do in fact roughly cancel in the marginal 
distributions. 

In §2 of this paper, we present our quasar sample from the SDSS. In §3 we derive the 
RLF of quasars as a function of redshift and optical luminosity. We examine the effects of 
K corrections and sample incompleteness in §4, and we give the discussion and summary in 
§5 and §6, respectively. Throughout the paper we use a A-dominated fla t cosmology with 
70 km s~^ Mpc-\ Qm = 0.3, and Qa = 0.7 (e.g. ISpergel et allbood ). 



2. The SDSS quasar sample 



The SDSS (lYork et al.l l2000l ) is an imaging and spec troscopic survey of the sky using 
a dedicated wide-field 2.5m telescope ( IGunn et al.ll2006l ). The ima ging is carried out in 
five broad bands , ugriz, spanning the range from 3000 to 10,000 A flFukugita et al.lll996l: 



Gunn et al.lll998r) . From the resulting catalogs of objects, quasar candidates (IRichards et al. 



2OO2I ) are selected for spectroscopic follow-up. Spectroscopy is performed using a pair of 
double spectrographs with coverage from 3800 to 9200 A, and a resolution A/ A A of roughly 
2000. The SDSS quasar survey spectroscopically targets quasars with i < 19.1 at low redshift 
{z ^ 3) and i < 20.2 at high redshift {z ^ 3). The low-redshift selection is performed 
in ugri color space, and the high-redshift selection is performed in griz color space. In 
addition to the optical selection, a SDSS object is also considered to be a primary quasar 
candidate if it is an optical point source located within 2'.'0 of a FIRST radio source. All 
SDSS magn itudes mentioned in this paper have b een corrected for Galactic extinction using 
the maps of ISchlegel. Finkbeiner. fc David (119981 ). 



The sample we used is from the SDSS Data Release Three (DR3: lAbazajian et al.ll2005l ). 
The quasar catalog of the pR3 consists of 46,420 objects with luminosities larger than 
Mi = —22 (ISchneider et al.l 120051 ). The area covered by the catalog is about 3732 deg^. 
We reject 4683 objects that are not covered by the FIRST survey, and we only use the 
quasars which were selected on their optical colors (i.e., the quasars with one or more of the 
following target selection flags: QSO_HIZ, QSO_CAP and QSO_SKIRT; see Richards et al. 
2002) to avoid the bias introduced by the FIRST radio selection. The final sample consists 
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of 31,835 optically-selected quasars from the SDSS DR3 catalog, and covers a redsliift range 
of < z ^ 5 and a luminosity range of —30 ^ A'Li < —22. 

To include both core-dominated (hereafter FRl) and lobe-dorninate d (hereafter FR2) 
quasars, we match our sample to the FIRST catalog (jWhite et al.l 119971 ) with a matching 
radius 30". For the quasars that have only one radio source within 30", we match them again 
to the FIRST catalog within 5" and classify the matched ones as FRl quasars. The quasars 
that have multiple entries within 30" are classified as FR2 quasars. The sample contains 
2566 FIRST-detected quasars, including 1944 FRl quasars and 622 FR2 quasars. 

We use the integrated flux density (/«„*) in the FIRST catalog to describe the 20 cm 
radio emission. The total radio flux density of each FR2 quasar is determined using all 
of the radio components within 30". We note that we have excluded those FR2 quasars 
whose separations between lobes are greater than 1'. In fact, FR2 quasars represent a small 
fraction of the SDSS DR3 catalog, and FR 2 quasars with diameters greater than 1' are 
even rarer (jde Vries. Becker, fc Whitd |2006| ) . These numbers are too small t o affect the 
statistics below. Therefore we do not use more sophisticated procedures (e.g. Ilvezic et al. 
2002[lde Vries. Becker, fc Whitelbooeh to select FR2 quasars. 



3. RLF of quasars as a function of redshift and optical luminosity 



We define a radio-loud quasar based on its R paramete r, the rest-frame ra tio of the flux 
density at 6 cm (5 GHz) to the flux density at 2500 A (e.g. IStocke et al.lll992l ). For a given 
quasar, we calculate its observed flux der isity /firm at rest-fr ame 6 cm from fi^t (if detected) 
assuming a power-law slope of —0.5 (e.g. Ilvezic et al.ll2004bl ): and we determine its observed 
flux density Ur^m at rest-frame 2500 A by fitting a model spectruna to the SDSS broadband 



photometry (IFan et al.ll200ll : iJiang et al.ll2006l : [Richards et al.ll2006l ). The model spectrum is 
a power-law continuu m (f„ = Au"") plus a serie s of emission lines extracted from the quasar 
composite spectrum (jVanden Berk et al.ll200ll ). We integrate the model spectrum over the 
redshifted SDSS bandpasses to compare with the observed magnitudes. The parameters a 
and A are determined by minimizing the differences between the model spectrum magnitudes 



m 



model 



and the SDSS photometry m' 



obs . 



X 



jj^model jj^obs 



cr. 



obs 



where 0"°''* is the estimated SDSS photometry error in the i^^ SDSS filter. We constrain a to 
be in the range —1.1 < a < 0.1, and only use the bands that are not dominated by Lyman 
forest absorption systems. Finally /2500 is computed from the power-law continuum using 
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the best-fit values of a and A, and tlie radio loudness R is obtained by 



-R — /ecm/ /: 



2500- 



(2) 



The absolute magnitude M2500 at rest-frame 2500 A is calculated from /; 



2500- 



The FIRST survey has a 5cr peak flux density limit of about 1.0 mJy (IBecker. White, fc Helfand 



19951 ). although this limit is not perfectly uniform across the sky. For a quasar detected by 
FIRST, we determine the relevant limit directly from the FIRST catalog, while for a quasar 
undetected by FIRST, we measure the limit at the position of the nearest radio source (usu- 
ally within 10'). We find that the medi an value of the limits is 0.98 mJy, which has alread y 
included the effect of "CLEAN bias" JSecker. White, fc Helfan"dl Il995l : Iwhite et al.lll997h . 
Only ~ 4% of the quasars have limits above 1.1 mJy, so we use 1.1 mJy as the FIRST 
detection limit for our sample. 

Many sources in the FIRST images are resolved. The resolution effect causes FIRST to 
beconie more incomplete for extended objects near the detection limit (IBecker. White, fc Helfand 
I995I : IWhite et al.lll997l ). Furthermore, FR2 quasars are more incomplete than FRl quasars 
for integrated flux densities. For example, a double-lobe radio source with two identical 
components suffers from incompleteness twice as high as a single-component source of the 
same total flux density. Figure 1 (provided by R. L. White, private communication) shows 
the FIRST completeness as a function of integrated flux density. The completeness is com- 
puted using the observed size distribution and rms values of integrated flux densities from 
the FIRST survey for SDSS quasars, and has included all effects mentioned above. Quasars 
with fint > 5 mJy have a completeness fraction gcomp ~ 1 (100%); while for a quasar with 
1.1 mJy < fint < 5 mJy, its gcomp is measured from the curve. To correct for sample in- 
completeness, we use the weight of I /gcomp when we calculate the numbers of radio- loud 
quasars. 



When quasars with i? > 10 are defined as radio-loud (e.g. iKellermann et al.l Il989l ). 
FIRST is able to detect radio-loud quasars down to i ~ 18.9 based on Equation [2l the K 
corrections we applied and the FIRST detection limit of 1.1 mJy. The left panel of Figure 2 
shows the redshift and absolute magnitude distribution of our sample. 

In flux-limited surveys, redshift and luminosity are artificially correlated, making it 
difficult to separate the dependence of the RLE on redshift or luminosity. To break this 
degeneracy, we divide the M2500-Z plane into small grids. RLFs in individual M2500-Z grids 
are calculated and presented as squares in the right panel of Figure 2, where the square for 
each subsample is located at the median values of M2500 and z in that subsample. The RLE 
declines with increasing redshift and decreasing luminosity. One can see the trend more 
clearly in Figure 3, in which we plot the RLE in three small redshift ranges and three small 
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log ( yZTrlp ) = feo + log(l + z)+ bM{M25oo + 26), (3) 



magnitude ranges. 

We assume a simple relation to model the RLF as a function of redshift and absolute 
magnitude, 

RLF \ 
7 

where b^, and 6a/ are constants. We use the RLFs calculated from the grids that include 
more than one radio-loud quasar, and use median values of z and M2500 in each grid. Sta- 
tistical uncertainties are estimated from Poisson statistics. The best fitting results found by 
regression fit are, 60 = -0.132 ± 0.116, b, = -2.052 ± 0.261, and bM = -0.183 ± 0.025. This 
implies that when RLF <^ 1, RLF oc (1 + z)~^'°^^L°p^^^, where Lopt is the optical luminosity. 
The of this fit is 52.8 for 47 degrees of freedom (DoF), and the confidence levels of bz 
and 6m are shown in Figure 4. The null hypothesis that b^ = and 6m = is rejected at 
> 5cr significance. The results are projected onto two-dimensional plots in the left panels of 
Figure 5, where filled circles are the RLFs calculated from individual grids in Figure 2 and 
dashed lines are the best fits. The upper panel shows the RLF as a function of redshift after 
correcting for the luminosity dependence. At M2500 = —26, the RLF drops from ~24.3% 
to ~4.1% as the redshift increases from 0.5 to 3. The lower panel shows the RLF as a 
function of luminosity after correcting for the redshift dependence. At z = 0.5, the RLF 
decreases rapidly from ~24.3% to ~5.6% as the luminosity decreases from M2500 = —26 to 
^2500 = —22. Therefore the RLF of quasars is a strong function of both redshift and optical 
luminosity. 

To probe whether the trend seen in Equation [3] is related to the radio-loud criterion 
adopted, we define a radio- loud quasar if R is greater than 30 instead of 10. In this case 
FIRST is able to detect radio-loud quasars down to i ~ 20.0. We calculate RLFs for the 
quasars with i < 20.0 using the same method illustrated in Figure 2, and model the RLF 
with Equation [3l The best fitting parameters are given in Table 1 and the confidence levels 
of bz and 6m are shown in Figure 4. Figure 5 shows the RLF as a function of redshift 
and luminosity. The relation gives similar results for both radio-loud criteria. We note 
that the sample of i < 20.0 at 2; < 3 is not complete since some bins in M25oo^z space are 
not sampled by SDSS. However, this incompleteness does not bias our results because we 



t (e.g. 


Peacock, Mi 


ler, & Loneair 


1986; 


Goldschmidt et al. 


19991). As we do for 



the i?-based RLF, we use two criteria to define radio-loud quasars: Lr (luminosity density 
at rest-frame 6 cm) > 10'^^ ergs s~^ Hz~^ and Lr > 10^^'^ ergs s^^ Hz~^. In the two cases, 
FIRST is able to detect radio-loud quasars up to 2; ~ 2.1 and 3.5, respectively. We model 
the RLF using Equation [3] and repeat the analysis. The best fitting results are shown in 
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Table 1, and Figures 4 and 6. The RLF based on Lr is correlated with z and M2500 in the 
same manner as the i?-based criteria. 



4. Effects of the K corrections and sample incompleteness 

When applying the K corrections, we assumed that the slope of the radio continuum is 
—0.5 and used a model spectrum to determine the optical continuum slope. To investigate 
the effect of the K corrections, we performed several experiments. First, for a given quasar 
at z, we calculated its /2500 and M2500 from the magnitude in the SDSS band whose effective 
wavelength is closest to 2500(1 + z) A assuming a slope of —0.5 (Test 1). As before, we 
corrected for the contribution from emission lines. Because the SDSS ugriz photometry 
covers a wavelength range of 3000 to 10,000 A, the K corrections for z < 3 require no 
extrapolation. Second, we assumed two extreme cases for the radio and optical slopes: in 
Test 2, we took the optical slope to be 0.0 and the radio slope as —1.0, while in Test 3, the 
optical slope was —1.0 and the radio slope was 0.0. The results of the fit to Equation [3] are 
listed in Table 1. The values of bz and Bm recalculated under these tests differ by less than 2a 
from the original values, and the null hypothesis that bz = and Bm = is rejected at > 5a 
significance in all these tests. Therefore the effect of the K corrections on our conclusions is 
small. 

We investigated the reliability of the relation described by Equation [3] for different 
definitions of radio loudness and for luminosities in different optical bands. For example, 
we defined the R param eter as the rest-frame ra tio of the flux density at 6 cm to the flux 



density at 4400 A (e.g. iKellermann et al.l Il989l ) instead of 2500 A, and we repeated the 



analysis in §3 (Test 4). In Test 5, we determined the RLF as a function of z and Mj (instead 
of M25oo)- Mi is the absolute magnitude in the rest-frame i band, and was calculated using 
the method described in §3. The results are listed in Table 1. In these cases the RLF is still 
strongly dependent on redshift and luminosity, and the relation described by Equation [3] is 
not sensitive to the details of how radio loudness is defined. 

When examining the dependence of the RLF on M2500 and z, we note that contours of 
constant RLF in Figure 2 roughly coincide with contours of constant apparent magnitude. 
This is illustrated in Figure 7(a), which shows that for different redshift bins, the relation 
between the RLF and apparent i magnitude is independent of redshift at i > 17.5. The RLF 
does decrease with redshift at z < 17.5, although with large error bars. This "conspiracy" 
of strong dependence of the RLF on apparent magnitude raises the concern that our results 
have been affected by flux- dependent selection effects. In this paper we are considering 
the RLF for optically-selected quasars, and the SDSS color selection is highly complete at 
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z < 2.2 flRichards et al.ll2006l ). so optical selection effects are not likely to seriously affect the 
RLF determination. However, the SDSS quasar sele ction becomes increasingly incomplet e 
for objects with very red int rinsic colors ( a < — 1.5, iFan et al.ll200ll : [Richards et al.ll2002l ). 
especially at high redshift. IWhite et al.l (120061 ) found a strong correlation between radio 
loudness and optical color using the SDSS sample. We reproduce this dependence in the 
upper panel of Figure 8. The RLF rises with increasing A{g — i), which is the difference 
between the observed g — i and the median g — i color of quasars at that redshift, following 
Hopkins et al.l (120041 ). To examine whether this RLF-color dependence affects the relation 
in Equation [3|, we divide the quasar sample into several A{g — i) bins, and calculate the RLF 
as a function of M2500 and z for each bin of intrinsic quasar colors. We find that although 
the average RLF increases toward redder continuum, the RLF is still a strong function of 
M2500 and z within each bin, similar to the relation shown in Equation [31 The lower panel 
of Figure 8 gives an example for the bin of —0.1 < A{g — i) < 0.1. Note that more than 
80% of the quasars in the sample lie in the range of— 0.3 < A{g — i) < 0.2, within which the 
RLF-color relation is relatively flat. Therefore, although we can not determine accurately the 
RLF evolution for the reddest few percent of quasars where SDSS is incomplete, the strong 
correlation between the RLF and both redshift and luminosity is not strongly affected by the 
RLF-color relation over the color range in which the SDSS selection is essentially complete. 

In order to examine radio selection effects, we performed the following tests. 



Did we miss FRl quasars due to the 5" radio catalog matching? We used a 10" matching 
instead, and found that the number of FRl quasars increases by only 3.7%. We also 
found that these additional sources increase the RLF by similar factors at both high 
and low redshift and both high and low luminosity, and thus have little effect on bz or 
bn- 

Did we measure the radio fluxes of F R2 quasars correctly? We compared FIRST with 
the NRAO VLA Sky Survey (NVSS; ICondon et al.lll998l ). which has a resolution of 
45". The FIRST and NVSS fluxes of most FR2 quasars are in good agreement. FIRST 
has a resolution of 5" and may overresolve radio sources larger than about 10". These 
sources are rare and very bright in radio (usually > 50 mJy), well above the radio- loud 
division in our analysis. 

Were there quasars detected by NVSS but not by FIRST? We matched our sample 
to the NVSS catalog within 15". We found that about 6% of the matched quasars 
were not detected by FIRST, and 80% of these additional sources are FRl sources 
with offsets more than 5" from the SDSS positions. They increase the RLF by similar 
factors at both high and low redshift and both high and low luminosity, and thus do 
not significantly change the trend in the RLF. 
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4. How did the incompleteness of FIRST at the detection hmit affect our results? In 
Section 3 we use 1.1 mJy as the FIRST detection limit, and we correct for the in- 
completeness caused by the resolution effect. Here we set two tests to examine the 
incompleteness near the FIRST limit. In Test 6 we use a limit of 1.5 mJy (the corre- 
sponding i ~ 18.6) to determine the RLF for R > 10; while in Test 7 we use a limit 
of 3 mJy (the corresponding i ~ 18.9) to calculate the RLF for R > 30. Note that the 
FIRST completeness measured from the completeness curve in Figure 1 is 75% at 1.5 
mJy and 95% at 3 mJy. The results of the tests are given in Table 1, and they show 
the similar trend of the RLF on redshift and optical luminosity as we obtained above. 

Therefore, based on these tests we conclude that the strong dependence of the RLF on M2500 
and z is highly likely to be physical, and the radio incompleteness and selection effects we 
have considered are not responsible for the dependence, though we cannot rule out other 
unexplored selection effects. 



Discussion 



Most of the previous studies of the RLF are based on samples of tens to hundreds of 
quasars. Considering that the RLF is only ~ 10% on average, these samples are not large 
enough to study the two-dimensional distribution of the RLF on M2500 and z. This makes 
it difficult to uncover the relation of Equation [3] using only the marginal distribution of the 
RLF. We calculate the marginal RLF as a function of M2500 and z for our sample, which 
is shown in Figures 7(b) and 7(c). One can see the marginal RLF is roughly independent 
of both M2500 and z, because the dependence on M2500 and z roughly cancel out due to the 
Mo^no-z degeneracy. This result is in quantitative agreement with the marginal RLF derived 
bv llvezic et aP J2002h . 



It has been suggested that high-redshift quasars show little difference in their rest-frame 
UV/optical and X-ray properties from those of low-redshift quasars. Their emission-line 
strengths and UV continuum shapes are very similar to those of low-redshift quasars (e.g. 
Barth et al.l l2003l : iPentericci et al.ll2003l : iFan et al.ll2004l ). the emission line ratios indicate 
solar or supersolar metallicity in emission-line regions as found in low-redshift quasars (e.g. 
Hamann fc Ferlandlll999l : iDietrich et al.ll2003l : iMaiolino et al.ll2003l ) , and the o ptical-to- X-ray 



2OO5I : ISteffen et al.l l2006l : IShemmer et al.l l2006l ). These measurements suggest that most 



3 opt 

flux ratios and X-ray cont i nuum shapes show litt le evolution with redshift (e.g. lStrateva et al. 



quasar properties are not sensitive to the cosmic age. However, Figure 2 shows that the RLF 
evolves strongly with redshift, thus the evolution of the RLF places important constraints 
on models of quasar evolution and the radio emission mechanism. 
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Equation [3] implies a strong correlation between t he RLF and opti cal luminosity. Using 
a large sample of low-redshift (0.03 < 2; < 0.3) AGNs, iBest et al.l (120051 ) find that radio- loud 
AGNs tend to reside in old, massive galaxies, and that the fraction of radio-loud AGNs 
is a strong function of stellar mass or central black hole mass (e.g., the fraction increases 
from zero at a stellar mass of 10^°Mq to 30% at a stellar mass of 5 x IO^^M^d). Assumin g 
that optical luminosity is roughly proportional to black hole mass ( jPeterson et al.ll2004l ). 
their radio-loud fraction of AGNs is also a strong function of optical luminosity. This is in 
quahtative agreement with Equation El although our quasars are more luminous than their 
low-redshift AGNs. 

Equation [3] also shows that t he RLF is a strong f unction of redshift. By stacking FIRST 
images of SDSS-selected quasars, IWhite et al.l (120061 ) recently found that the median i? is a 
declining function with optical luminosity. After correcting for this effect, they claimed that 
the median R is independent on redshift, which seems inconsistent with our result that the 
RLF is a strong negative function of redshift. However, the median R and the RLF are not 
identical. The median R is determined by the majority of quasars with low R values (i.e. 
radio-quiet quasars); while the RLF is the fraction of quasars exceeding a threshold in R, and 
therefore corresponds to the behavior of the small fraction of quasars with high R values (i.e. 
radio-loud qu asars). There are two natural way s to interpret the evolution of the RLF with 
redshift (e.g. iPeacock. Miller, fc Longairl Il986l ): (1) This may be due to the cosmological 



evolution of quasar radio properties, such as R and Lr- For instance, a decreasing R results 
in a decreasing RLF for increasing redshift. (2) This could be simply caused by the density 
evolution of different popula tions of quasars (e.g. radio-loud and radio-quiet quasars). The 
results of IWhite et al.l (120061 ) may have ruled out the first explanation and leave the second 
one: the different density evolution behaviors for the two classes of quasars. For instance, 
there are more radio- loud quasars at low redshift, but the fraction of radio-loud quasars is 
small, so they do not change the median R, which is still dominated by radio-quiet quasars. 
This claim is based on stacked FIRST images. To distinguish between the two explanations, 
one needs to determine the radio luminosity function of quasars in different redshift ranges, 
including the radio-quiet population, going to radio fluxes muc h fainter than those p robed 
by FIRST. Deep surveys such as the Cosmic Evolution Survey (jSchinnerer et al.ll2004l ) that 
cover a wide redshift range and reach low luminosity in both optical and radio wavelengths 
are needed to interpret the evolution of the RLF. 
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6. Summary 

In this paper we use a sample of more than 30,000 optically-selected quasars from SDSS 
to determine the RLF of quasars as a function of redshift and optical luminosity. The sample 
covers a large range of redshift and luminosity. We study the RLF using different criteria to 
define radio-loud quasars. After breaking the degeneracy between redshift and luminosity, 
we find that the RLF is a strong function of both redshift and optical luminosity: the RLF 
decreases rapidly with increasing redshift and decreasing luminosity. The relation can be 
described by a simple model, given by Equation [31 We have done a series of tests to examine 
the impact of flux-related selection effects, and find that the dependence of the RLF on 
redshift and luminosity is highly likely to be physical. 

The RLF is one of a few quasar properties that strongly evolve with redshift, so the evo- 
lution of the RLF places important constraints on models of quasar evolution and the radio 
emission mechanism. By comparing our results with the behavior of the median R derived 
from stacked FIRST images, we find that the evolution of the RLF with redshift could be 
explained by the different density evolution for radio-loud and radio-quiet quasars. Substan- 
tially deeper wide-angle radio surveys which obtain fluxes for the radio-quiet population are 
needed to fully understand the physical nature of this evolution. 
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Table 1. Results of fits to Equation [3] 
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-0.132 ± 0.116 


-2, 


.052 ± 0. 
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.025 
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Fig. 1. — The FIRST completeness as a function of integrated flux density (provided by R. 
L. White, private communication). The completeness is computed using the observed size 
distribution of SDSS quasars and rms values of integrated flux densities from the FIRST 
survey. 
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Fig. 2. — Left panel: Redshift and absolute magnitude distribution for 20,473 quasars with 
i < 18.9 in our sample. The M25oo^z plane is divided into small grids to break the redshift- 
luminosity dependence. RLFs of quasars are calculated in individual grids. Right panel: 
The i?-based RLFs in individual M2500-Z bins. The square for each subsample is positioned 
at the median values of M2500 and z in that subsample. The RLF of quasars declines with 
increasing redshift and decreasing luminosity. 
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Fig. 3. — RLF in three small redshift ranges and three small magnitude ranges. Dotted 
hnes are the best model fits. 
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Fig. 5. — The i?-based RLF as a function of z and M25oo- The upper panels show the RLE 
as a function of redshift after correcting for the luminosity dependence and the lower panels 
show the RLF as a function of luminosity after correcting for the redshift dependence. Filled 
circles are RLFs calculated from individual M2%m-z grids and dashed hnes are the best model 
fits. Poisson errors are also given in the figure. 
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Fig. 6. — The L^-based RLF as a function of z and M25oo- The upper panels show the RLE 
as a function of redshift after correcting for the luminosity dependence and the lower panels 
show the RLF as a function of luminosity after correcting for the redshift dependence. Filled 
circles are RLFs calculated from individual M2500-Z grids and dashed hnes are the best model 
fits. Poisson errors are also given in the figure. 
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Fig. 7. — (a): RLF as a function of apparent magnitude i for five redshift bins. The RLE 
declines with increasing and the curves for different redshift bins follow the same RLF-i 
dependence, (b) and (c): Marginal RLF as a function of redshift and luminosity for our 
sample. The marginal RLF is roughly independent of both redshift and luminosity due to 
the M250Q-Z degeneracy. 
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Fig. 8. — Upper panel: RLF as a function of A{g — i). The RLE rises with increasing 
A{g — i) at z < 2.2. Lower paneh RLF as a function of z and M2500 in a small color range 
of -0.1< A(c/-z) <0.1. 



